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Abstract 

We consider the spin-i anisotropic XY chain in a transverse (z) field with the Dzyaloshinskii-Moriya 
interaction directed along z-axis in spin space to examine the effect of the Dzyaloshinskii-Moriya interac- 
tion on the zz, xx and yy dynamic structure factors. Using the Jordan- Wigner fermionization approach 
we analytically calculate the dynamic transverse spin structure factor. It is governed by a two-fermion 
excitation continuum. We analyze the effect of the Dzyaloshinskii-Moriya interaction on the two-fermion 
excitation continuum. Other dynamic structure factors which are governed by many-fermion excitations 
are calculated numerically. We discuss how the Dzyaloshinskii-Moriya interaction manifests itself in the 
dynamic properties of the quantum spin chain at various fields and temperatures. 



1 Introduction. Jordan- Wigner representation 

Recently the dynamic properties of quasi-one-dimensional quantum spin systems have become an intense 
area of research. On the one hand, a variety of quasi-one-dimensional materials were discovered during 
last decades. On the other hand, many exactly solvable statistical mechanics models refer to one spatial 
dimension. Very often these models permit to make a detailed analysis not only of their thermodynamic 
properties but also of their dynamic properties. Such theoretical investigations are necessary for interpreta- 
tion of the experimental data observable in the scattering or resonance experiments on quasi-one-dimensional 
compounds pQ. 

In this paper we analyze the effect of the Dzyaloshinskii-Moriya interaction on the dynamic properties of 
quantum spin chains. The Dzyaloshinskii-Moriya interaction is present in many low-dimensional materials 
[2ll3l[ll[5l[6l[7l[8l[9l [TUl HU H2l H3J HH [15]. Such antisymmetric exchange interaction, D nm • [s„ x s m ], 
takes place, if allowed by crystal symmetry, due to spin-orbit coupling |16j and is weaker than the symmetric 
Heisenberg superexchange interaction J nm (s n • s m ). In spite of this, this interaction has a number of impor- 
tant consequences and may cause a number of unconventional phenomena. Interestingly, the Dzyaloshinskii- 
Moriya interaction may appear while analyzing the noncquilibrium steady states of quantum spin chains 
with currents [17\ \18\ I19j . It may also appear in the quantum spin representation in the stochastic kinetics 
of adsorption-desorption processes [20] (see also [2Pj). 
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The analysis of the effect of the Dzyaloshinskii-Moriya interaction on the dynamic quantities of quantum 
spin chains was reported in several papers [SUES]. In particular in Refs. [22j[23] such analysis was performed 
using the symmetry arguments and field-theoretical methods for the isotropic Heisenberg (XXX) chain. 

In our study we restrict ourselves to a simpler model, i.e. anisotropic XY chain, the dynamic properties 
of which are amenable to detailed analytical and numerical analysis. We should also note that CS2C0CI4 
provides a new example of spin-i XY chain [24] that may renew interest in the calculation of observable 
quantities for such a chain 25, 26, 27 . 

To be specific, we consider TV spins one-half governed by the following Hamiltonian 

H = E ( JXs n S n+l + J V <sl + l) +Y, D (< S «+1 " S « S «+l) + E ^ ( L1 ) 

n n n 

Here J x , J v are the anisotropic XY exchange interaction constants, D is the z-component of the Dzyaloshin- 
skii-Moriya interaction and f2 is the transverse magnetic field. Such model was introduced in Refs. [28l [29] , 
Some of its dynamic properties were examined in Refs. [301 Ell 1321 OS [34]. In particular, the transverse 
dynamic susceptibility Xzz(k,uj) of the model (jl.ip was derived explicitly for n = [30] and k ^ [32]. 
Moreover, the case of isotropic XY interaction J x = J y (in this case the Dzyaloshinskii-Moriya interaction 
can be eliminated from the Hamiltonian by a spin axes rotation) was analyzed in some detail [34] 135] . 
The performed analysis is based on the Jordan- Wigner transformation, 

d = si, c ri = (-2sf)(-24)...(-2<_ 1 )s-, n = 2,...,N, (1.2) 

which transforms (jl.ip into the following bilinear Fermi form 

J + iD , J —ID , 7/ + + x n ( + 1 



H = 2_j C n Cn+l ^ C " C «+1 + 2 ( C « C ™+1 ~ C « C ™+0 + Q [°n C n ~ 3 J J ( L3 ) 

with J = \ (J x + jv) and 7 = i (J* - J*). 

In our calculations we consider both periodic and open boundary conditions. Of course, in the thermo- 
dynamic limit the results for bulk characteristics should be insensitive to the boundary conditions implied. 
In the former case, i.e. when s^+i = s", the bilinear Fermi form (|1 .3|) is periodic or antiperiodic depending 
on whether the number of fermions is odd or even. After the Fourier transformation, 

\= ^2 ex P c n, c n — — != ^2 ex P (— c K (1.4) 



N ~ VN 

(here k takes N values ■jj-n for periodic boundary conditions or ^ (n + i) for antiperiodic boundary con- 
ditions; n = -f , -f + 1, . . . , f - 1 for N even and n = -^f 1 ,-^ + 1, • • ■ , ^f 1 for ^ odd )> and thc 
Bogolyubov transformation with Fermi-operators (3 

(3 K = iu K c K + v K c± K , c K = -iu K (3 K + v K f3t K (1.5) 

with 



. . 1 / O, + J COS K 1 / O + J COS K 

. K =sgn( 7 sm K )-^l + ^--— , « K = ^1 — , 



A K = a/ (f2 + Jcosk) 2 + 7 2 sin 2 k (1.6) 
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the Hamiltonian (|1.3p becomes 

Tl — , , , , , _ 



A K = L>sinK + A K . (1.8) 

It should be noted here that only the energy spectrum A K (|1.8p but not the coefficients of the Bogolyubov 
transformation u K , v K (|1.6p depends on £>. Using (|1.8|) one immediately finds that the energy spectrum is 
gapless when j 2 < D 2 and f2 2 < J 2 + D 2 — 7 2 or when j 2 > D 2 and f2 2 = J 2 . 

In our numerical calculations we use open boundary conditions 36J. The bilinear Fermi form (|1.3|) with 
open boundary conditions can be brought into the form (|1.7j) after a linear transformation 

JV N 

Vk = Y^ i.9k 3 Cj + h kj cf) , Cj ; = ^ (g^rik + h k jVk) ( L9 ) 



where 



j=i k=i 



N N 



^ (gkjAjm — hkjB* m ) — Akgkm, ^2 (gkjBjm — hkjA* m ) — Akhkm, 

iV iV 

(fffciffgj + h k jh* qj ) = Skq, ^2 (gkj9km + h* kj h km ) = d jm , 
j=i j'=i 

iV N 

(9kjh qj + ftfcjfirgj) = ^2 (gkjh* km + h* k3 g km ) = (1.10) 
i=i fc=i 

^<W + i (J + iD) $ m ,n+i + 5 (•/ - lD) ^m,n-i 5 Bnm = §7 (<*m,n+i - ^m.n-i)- We are not aware 
of a general analytical solution for this problem. For the two particular cases, namely, the anisotropic XY 
chain without field and the Ising chain in a transverse field such solutions can be found e.g. in Ref. [37] . In 
our study we solve Eqs. (|1.10p numerically for chains of about a few hundred sites. 

The relation between the spin model (jl.ip and the noninteracting Jordan- Wigner fcrmions (jl.7p is a key 
step in the statistical mechanics calculations for one-dimensional spin-i XY systems. 
In our study of the dynamic properties we focus on the dynamic spin structure factors 

S aa (K,u}) = — ^^%xp(i/tn) / dtcxp (iwt) ((s^(t)sf +n ) - (sf)(sf +n )) , a = x,y,z. (1.11) 

JV j=ln=l *'-«= 

The dynamic structure factors are of considerable importance since they are directly comparable with in- 
elastic neutron scattering experiments of some quasi-one-dimensional substances. The dynamic transverse 
spin structure factor S zz (k, lj) can be easily evaluated analytically (Section 2). The transverse dynamics is 
governed exclusively by a two-fermion excitation continuum the properties of which in the case of XY chain 
without the Dzyaloshinskii-Moriya interaction were discussed earlier |38] . Another quantity which is also 
governed by the two-fermion excitation continuum is the dynamic dimer structure factor [39, 40, 41J. There- 
fore, the effect of the Dzyaloshinskii-Moriya interaction on the two-fermion excitation continuum deserves a 
separate discussion (Section 3). The xx and yy dynamic structure factors are computed numerically (Section 
4) . We compare and contrast the properties of different dynamic structure factors at different values of the 
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transverse field and temperature emphasizing the effect of the Dzyaloshinskii-Moriya interactions. We end 
up with conclusions (Section 5). 

Some preliminary results of this study were announced in the conference paper [42] . 



2 Dynamics of transverse spin correlations 

We start with the zz dynamic structure factor of the model (|l.ip . For analytical calculation of this quantity 
one can consider only periodic (or only antiperiodic) boundary conditions for the bilinear Fermi form (|1.3j) . 
Using the relation between spin and Fermi operators (|1.2[) and the transformations (|1.4[) and (|1.5p . (| 1 . 6[) 
after standard calculations using the Wick-Bloch-de Dominicis theorem we arrive at 

S zz (k,u>) = / d/«i ((u 2 Kl u 2 Kl _ K - u Kl u Kl - K v Kl v Kl - K ) n K1 (1 - n Kl - K ) 5 (uj + A Kl - A Kl _ K ) 

J — 7T 

+ Kj^-k + u Kl u Kl - K v Kl v Kl - K ) n Kl n- Kl+K 5 (uj + A Kl + A_ Kl+K ) 
+ + u Rl u Kl - K v Kl v Kl - K ) (1 - U-kJ (1 - n Kl - K ) 5(uj- A_ Kl - A Kl _ K ) 

+ (vl x v 2 , 1 _ K - u Kl u Kl - K v Kl v Kl - K ) (1 - 7i_ Kl ) n- Kl+K S (ui - A_ Kl + A_ Kl+K )) (2.1) 

where n K = (exp(/3A K ) + 1) 1 is the Fermi function. This result agrees with the corresponding formula for 
the transverse dynamic susceptibility Xzz(^,^) derived in Ref. [32 . 
Introducing the function 

(Q -j- J cos ( Kl - I )) (O + Jcos ( Kl + | )) - 7 2 sin ( Kl - |) sin ( Kl + |) 

A «l"f A «l + f 

the dynamic structure factor S zz (k,uj) (|2.1[) can be expressed as follows 

S zz (k,uj) = J d«i ^ + ^ _ nrei _^) n« 1+ *<5 (w - A Kl _| +A Kl+ s) 

+ ((1 ~ n^-f) (1 ~ ) * (- - A K1 - S - A_ B1 _ f ) 

+n K1+f n_ K1 + - 5 (w + A K1+f + A_ K1+f ))) . (2.3) 

In the limit of isotropic XY interaction 7 = Eq. (|2.3p yields the result obtained earlier [34]. In the limit 
D = (when A K = A K = A_« > 0) and T = ((3 -> 00) Eq. ([23]) becomes 

S„(«,w) = I" d Kl 1 = IpA s (uj - A B1 _ a - A K1+f ) . (2.4) 

This coincides with the expression obtained earlier in Ref. [38] . We notice that Eq. (|2.4[) is more generally 
valid for the case 7 2 > D 2 (when A K > 0) and T = 0. In the case D 2 > j 2 and T = or in the most general 
case of arbitrary D and T > (as well as in the case D = but T > which was not considered in Ref. 38J) 
the zz dynamic structure factor exhibits new qualitative features in comparison with the analysis reported 
in Ref. [38] , Again the zz dynamic structure factor is governed exclusively by two-fermion excitations as 
can be seen from Eq. (12. 3|) , however, for D 2 > j 2 , T = or for T > all three (^-functions in Eq. (|2 .3|) may 
come into play. 

In Fig. 1 we show the gray-scale plots for S zz (k,u>) at low temperature (ft J = 50) for several typical 
sets of parameters (J = 1, 7 = 0.5, D = 0, 0.5, l,O = 0, 0.5, 1, ^-). For 7 2 > D 2 (right column in 
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Fig. 1) at such a low temperature only one two-fermion excitation continuum is relevant (the one which 
arises owing to 8 (oj — A Kl _* — A_ Kl _s) in (|2.3j) ) whereas in the opposite case D 2 > j 2 (left column in Fig. 
1) all three two-fermion excitation continua contribute to transverse dynamics. In Fig. 2 we demonstrate 
typical low-temperature frequency profiles of S zz {0,ui), S zz (^,uj), S zz (ir,uj) for a chain with J = 1, 7 = 0.5, 
= 0.5 and different values of the Dzyaloshinskii-Moriya interaction D = 0, 0.5, 1. In Fig. 3 we show 
the gray-scale plots for S zz (k,uj) at intermediate and high temperatures (/3 J = 10, 1, 0.1) for a chain with 
,7=1,7 = 0.5, D = 1,Q, = 0.5. All these data are used herein to discuss generic properties of a two-fermion 
dynamic structure factor (Section 3) and to compare different dynamic structure factors (Section 4). 

3 Two-fermion excitation continua 

The zz dynamic structure factor (|2.3[) can be rewritten in the form 

S zz (k,oj) = J2 S zJ{^^): S { zJ(k,oj)= f dmB^im.^jC^iKu^sL-E^iK^K)) (3.1) 

j=l J -' K 

with 

B^( Kl , K )=B^( Kl ,n)= 1 - f [ Kl ' K \ 

B( 2 \ Kl , K )= 1 -±Ip^, (3.2) 

= (1 - (1 - "-«!-$) . 

C (2) (k 17 k) = (l - n K1 _%) n K1+ §, 

C (3) (ki,k) = n Kl+ « n_« 1+ £ , (3.3) 



-Kl- 



£ (1) (ki,k) = A Kl _ f +A_ 
£ (2) (ki,k) = A Kl _ f - A Kl+f , 
£( 3 )( Kl , K ) = -A Kl+s -A_ rei+f . (3.4) 

In accordance with (|3.1[) we distinguish three two-fermion excitation continua (they correspond to j = 1, 2, 3 
in Eq. (|3.ip ) which govern S zz (k, cj). The gray-scale plots of S'li^rt, cj), j = 1,2,3 for a typical set of 
parameters, J = 1, 7 = 0.5, D = 1, f2 = 0.5, /3 = 50, are shown in Fig. 4. As can be seen from (|3.1| and 
(13. 2|) , (|3.3p , (|3.4p the specific features of the considered two-fermion dynamic structure factor which describes 
the dynamics of the transverse spin fluctuations are controlled by the S^-functions. Therefore, in Fig. 5 
we also plot S^z (k, u>) ()3.1|) for the same set of parameters as in Fig. 4, although, with B^\ki,k) = 1. 
By comparing Figs. 4 and 5 one can distinguish the specific features (coming from the -functions J3.2[l ) 
and the generic features (coming from the C^-functions (|3.3p and the i?W-functions (|3.4p ) of the dynamic 
quantity considered (see discussion below). 

We remark that although we were not able to find a simple analytical form of the most important lines in 
the k—u> plane characterizing the two-fermion excitation continua for a general case of the spin chain (jl.ip it 
is easy to determine these functions numerically for any set of parameters using MAPLE or/and FORTRAN 
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codes. We also recall that in the simplest case of the isotropic XY model in a transverse field (7 = 0, 
D = 0) the corresponding results have been derived analytically [43]. However, in the case 7 7^ 0, D = 
the analytical results have been reported only in the limiting cases Q = or 7 = 1 [38]. Of course, the case 
7 =/= 0, D =/= considered in the present paper is even more complicated. In what follows we take a typical 
set of parameters J = 1, 7 = 0.5, D = 1, f2 = 0.5. 

We begin with the case of infinite temperature T — ► 00 (/3 = 0). The two-fermion dynamic structure 
factor may have nonzero values in the plane wave- vector k - frequency u> if the equation 

lu - £; (j) (k!,k) = (3.5) 

has at least one solution k*, —it < k* < it. In Fig. 6 (panels a (j = 1), b (j = 2), c (j = 3)) we show 
the regions in the k-lu plane in which equation (|3.5p has four solutions (dark-gray regions), two solutions 
(gray regions) or has no solutions (white regions). The bounding lines of the regions in which equation (|3.5[) 
has solutions constitute the upper (uu\k)) and the lower (lu^\k) < lju\k)) boundaries of the two-fermion 
excitation continuum, respectively; moreover, for some regions of the wave-vector k the lower boundaries 
u>i (k) may be equal to zero. Alternatively, we may find the upper and the lower boundaries of the two- 
fermion excitation continuum seeking for the maximal and minimal values of (k\, k) while n\ varies from 
— 7r to 7r, i.e., 



(k) = max (eW(ki,k)\, U j) (K) = min (0, £ w (/ci, k)\ . (3.6) 



We have found that u)v\k) and oj^\k) ^ occur for the values of «i, K*, which satisfy the equation 



OKI 



0. (3.7) 



Moreover, Eq. (|3.7p also holds along the boundary u>(k) between the gray and the dark-gray regions in the 
upper (j — 1) and the middle (j = 2) panels in the left column in Fig. 6. On the other hand, the quantities 



S U) («, oj) = f dmS^ (ki , k)S (oj - (ki, «)) = 



(3.8) 



where {k*} are solutions of Eq. (|3.5[) , may exhibit a van Hove singularity along the line ojg (k) — E^ (k*, k), 
where n* satisfies Eq. Q3.7p . Thus, the mentioned boundary lines in the panels in the left column in Fig. 6 
are the lines of van Hove singularities akin to the density of states effect in one dimension. Further, we have 
found that for almost all cases -r^-zE^ (ki, k) ^ for the values of Ki which satisfy (|3.5[) with oj — u>i' (k) 
that obviously implies a familiar square-root divergence 

SW(K,u)oce-$ (3.9) 

when oj approaches w, (k), e = u> — u>s (k) ■ However, for j — 2 and k ps 1.07844531 (see panel b in 
Fig. 6) Eq. ([37]) holds for m « 2.16480069 and at this point -£^E^ ( Klj K ) = but J^eW(ki,k) » 
2.96548741 ^ 0. As a result, if oj is in the e- vicinity of u4 2) (k w 1.07844531) w 0.78594502 the quantity 
S( 2 \k,uj) (|3.8[) exhibits singularity with another exponent 

S (2) (/t re 1.07844531, 0.78594502 ± e) oc e"i. (3.10) 
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We illustrate different types of singularities in Fig. 7. In particular, in Fig. 7a one can see the square-root 
divergencies Q3.9j) , whereas in Fig. 7b aside from the square-root divergences Q3.9[) (solid and dotted curves) 
one can also see the dependence (|3.10l) (dashed curve). We notice that the e~i singularity remains for other 
values of D and is also present when D — > 0. For J = 1, 7 = 0.5, D = 0, O = 0.5 it occurs at k *=> 1.68213734 
while lo approaches ui s (k rss 1.68213734) = 0. Interestingly, this fact could not be detected in the earlier 
study on the zz dynamics in the anisotropic XY chain in a transverse field without the Dzyaloshinskii-Moriya 
interaction [38] since that study refers to the zero-temperature case when only the continuum j — 1 manifests 
itself in the transverse spin dynamics (see discussion after Eq. (|2.4p ). We also notice that the observation of 
the e~i singularity may be difficult because of the fact that this peculiarity takes place only at one specific 
value of k (in contrast to the e _ 2 singularity). However, for the values of wave- vector in the vicinity of this 
specific value one easily distinguishes a reminiscence of the e~3 singularity (see the dotted curve in Fig. 7b). 
Finally we note that for some of the lines characterizing the two-fermion excitation continua we can give 
simple analytical expressions. Thus, for j = 1 the maximum/minimum occurs at K\ = and 

K± = —it and hence the corresponding boundary lines are given by E<U(0,k) and —it, k). We did not 
find simple analytical expressions for the boundary lines between the white and the dark-gray regions and 
for the nonzero lower boundary (see panel a in Fig. 6). For j = 3 the upper boundary is given by E^{0,k) 
and £( 3 )(-7T,k). 

Next we turn to the zero-temperature case T = ((3 — > 00) for which the effect of the Fermi functions 
involved in C^)- functions becomes important. The Fermi functions contract a region of possible values of 
k,\ in Eq. (|3 . 5[) : now k% varies only within a part of the region between — tt and tt where C^\k\, k) 7^ 0. In 
Fig. 6 (panels d (j = 1), e (j = 2), f (j = 3)) we show the regions in the k-lu plane in which Eq. (|3.5[) taking 
into account the condition C^\k\,k) 7^ has four solutions (dark-gray regions in panel d), two solutions 
(gray regions), one solution (light-gray regions in panel e) or has no solutions (white regions). The upper 
and the lower boundaries may be also found according to Eq. (|3.6p . although with Ki varying within a part 
of the region between — tt and tt where C&(ki, k) 7^ 0. The values of k at which Wj (k) = correspond 
to the potential soft modes k c (see the right panels in Fig. 6). (Note that the soft modes may not occur 
owing to B^ -functions (compare Figs. 4 and 5).) The two-fermion dynamic quantities may exhibit the 
above discussed van Hove singularities if their occurrence is not prohibited by the Fermi functions (and 
B^' -functions, compare Figs. 4 and 5). Moreover, along the bounding lines between the regions of K—u> 
plane corresponding to different numbers of solution of Eq. (|3.5[) the two-fermion dynamic quantities may 
abruptly change their values exhibiting a finite jump (for example, we refer to the boundary line between 
the gray and light-gray regions in Fig. 6e; see also the panels b in Figs. 4 and 5). We can write down the 
analytical expressions for some characteristic lines seen in Figs. 6d, 6e, 6f. We introduce 



cosk± = ¥-i T9 ' 3 '- '-, flsin K± <0; 3.11 

J + D z — 7 Z 

K± obey the equation A K± = 0. For the chosen set of parameters J = 1, 7 = 0.5, D = 1, Cl = 0.5 we have 
k + w —1.24466864, w —2.67211739. The lines A_ K± _ K form the lower boundary of the continuum j = 1; 
the lines A K± _ K , — A K± + K form the lower boundary of the continuum j = 2; the lines — A re _ re± form the lower 
boundary of the continuum j — 3. The soft modes may occur at the values of wave- vector — 2«±, — k+ — K— 
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(j = 1), 0, ± (k+ - k_) - 2), 2/t±, k+ + K _ (j = 3). 

Finally, we emphasize a role of £>( J ) -functions (13. 2| which are responsible for the specific features of 
the dynamic transverse spin structure factor S zz (k,uj). The functions B^'(kx,k) modify and add some 
additional structure to S zz (n,tu) in the K—u) plane (compare Fig. 4 and Figs. 5, 6d, 6e, 6f referring to the 
low-temperature limit and Fig. 3c and Figs. 6a, 6b, 6c referring to the high-tcmpcraturc limit). In particular, 
the function B^(ki, k) removes the soft modes at k = ± (k + — k_) but not at K = from S zz (k,uj) (see 
Fig. 4b). Furthermore, comparing Figs. 4a, 4c and 5a, 5c one sees that van Hove singularities along the 
lines uj = E^(0, k), ui = E^(—n, «) (panels a) and along the lines uj — E^(0, «), ui = E^{—it, k) (panels 
c) disappear since B^ (0, n) = B^ (0, k) = (-tt, k) = B (3) (-tt, k) = 0. 

To summarize this Section, the two-fermion dynamic structure factors have a nonzero value only in a 
restricted area of the k—ui plane (two-fermion excitation continua) and may exhibit the van Hove singularities 
not only with exponent | but also with exponent |. Moreover, at zero temperature the two-fermion dynamic 
structure factors may exhibit jumps at which their values abruptly increase by a finite value. 

4 xx and yy dynamic structure factors 

In the present Section we calculate the xx and yy dynamic structure factors. For numerical calculations it 
is convenient to rewrite p. lip in the following form 



/oo 
dtexp(i^)((,s«W Sj Q )-(.s tt ) 2 ) 
-00 




(we have omitted jjJ2f=i m CHU) an d have used the relation (s"(— t)s" +n ) = (s"(t)s"_ n )*). As a re- 
sult, to get S aa (n,oj) according to (|4.ip we have to compute the time-dependent spin correlation functions 
(Sj(t)sj +n ), n = 0, 1,2, . . . with t varying from — oo to oo. We compute (s"(t)s" +n ) a = x,y numerically. 
For this purpose we express the spin operators entering the time-dependent correlation functions in terms of 
auxiliary operators tpr = c ; + ± c/ which according to Eq. (|1.9[) are linear combinations of operators rf^ , r/k- 
Since the x and y spin components at each site are essentially nonlocal objects in terms of Jordan- Wigner 
fermions (see (|1.2p ) the resulting expressions for the time-dependent spin correlation functions (s"(i)s" +n ), 
a = x,y are complicated averages of products of a large number of Fermi operators attached not only the 
sites j and j + n but to two strings of sites extending to the boundary of the chain. After applying the 
Wick-Bloch-de Dominicis theorem one gets an intricate result which can be compactly written as the Pfaf- 
fian of the 2(2j + n — 1) x 2(2j + n — 1) antisymmetric matrix constructed from elementary contractions 
(<Pl (t)(p^ n ), (ipf (t)<fm) i (^r WVm)) (f7 (%m) ■ The elementary contractions are easily expressed in terms 
of Afe, gkj, hkj (|1.9p . (jl.lOp . Finally we numerically evaluate the Pfaffians. Typically we take TV = 400, 
assume j = 31, 71, 131 and calculate the correlation functions with n up to 50, 70, 140 in the time range 
-50 50, -100 100, -300 300, -600 -r- 600. To remove the effect of a finite time cut-off we multi- 
ply the integrands in (|4.ip by exp(— ie\t\) with e = 0.01, 0.02. To be sure that our results pertain to the 
thermodynamic limit we examine in detail different types of finite size effects (for further details sec Rcfs. 
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In Fig. 8 (Fig. 9) we plot S xx (k,uj) (S yy (n,uj)) at low temperature (j3J = 50) for several typical sets 
of parameters (J = 1, 7 = 0.5, D = 0, 0.5, 1,0 = 0, 0.25, 0.5, 1). (For larger values of the transverse 
field SI there are no qualitative changes in these dynamic structure factors: they have a simple single-mode 
structure.) In Fig. 10 (Fig. 11) we show typical low-temperature frequency profiles of S xx (k,uj) (S yy (n,uj)) 
at K = 0, K = ^, K = 7r for a chain with J = 1, 7 = 0.5, SI = 0.5 and different values of the Dzyaloshinskii- 
Moriya interaction D = 0, 0.5, 1. In Fig. 12 (Fig. 13) we plot S xx (k,uj) (S vv (k,,ui)) at intermediate and 
high temperatures (/3J = 10, 1, 0.1) for a chain with J = 1, 7 = 0.5, D = 1, SI = 0.5. 

In contrast to transverse dynamic structure factor, the xx and yy dynamic structure factors are essen- 
tially more complicated quantities within the Jordan- Wigner method. Really, owing to a nonlocal relation 
between the x and y spin components and Fermi operators (|1 .2|) the xx and yy time-dependent spin corre- 
lation functions are expressed through many-particle correlation functions of noninteracting Jordan- Wigner 
fermions. Let us now discuss the obtained numerical results. First of all we note that both dynamic structure 
factors S xx (k,u>) and S yy (n,uj) show similar behavior for the taken value of 7 = 0.5; obviously they become 
identical in the isotropic limit 7 = 0. We start with the dynamic structure factors at low temperatures. 
As can be seen in Figs. 8 and 9 (and Figs. 12a and 13a) these dynamic structure factors show several 
washed-out excitation branches which are roughly in correspondence with the characteristic lines of the 
two-fermion excitation continua (compare three dynamic structure factors in the panel a and the panels b 
and c in Fig. 14; note that these quantities are shown for k that varies from — tt to 37r). Thus, although 
S xx (k,u>) and S yy (K,uj) are many-particle quantities within the Jordan- Wigner picture and hence they are 
not restricted to some region in the k—u> plane, their values outside the two-fermion continua are rather 
small. This observation (i.e. two-particle features dominate many-particle dynamic quantities S xx (k,u>) and 
S yy (n,uj) at low temperatures) agrees with our previous studies on isotropic XY chains 36, 34 (see also 
Ref. [38]). The constant frequency scans for several values of the wave- vector displayed in Figs. 10, 11 
show the redistribution of spectral weight S xx (k,u>) and S yy (n,u>) as the Dzyaloshinskii-Moriya interaction 
D increases. We note that these frequency profiles exhibit one or several peaks that may be relatively sharp 
or broad. The Dzyaloshinskii-Moriya interaction affects the positions of the peaks, their shapes and even 
their number (see, for example, the dependences S xx (0, to) vs to and S yy (0, to) vs u> displayed in Figs. 10a and 
11a). Constant frequency (and wave- vector) scans can be obtained for quasi-one-dimcnsional compounds by 
neutron scattering or resonance techniques and our findings may be useful in explaining the experimental 
data for the corresponding materials. As can be seen from our results, the Dzyaloshinskii-Moriya interaction 
clearly manifests itself in the frequency or wave-vector profiles of the dynamic structure factors that can be 
used in determining the magnitude of this interaction. 

It should be remarked that the dynamic structure factors of quantum spin chains are often examined 
within the framework of a bosonization approach [44, 23J. Note, however, that field-theoretical approaches do 
not apply to small length scales and short time scales when the discreteness of the lattice becomes important. 
Therefore, since these methods can describe only the low-energy physics, the high-frequency features nicely 
seen in Figs. 8 and 9 cannot be reproduced by these theories. 

As temperature increases, the low-temperature structure gradually disappears and the dynamic structure 
factors S xx (k,lo) and S yy (n,uj) become ^-independent in the high-temperature limit (see Figs. 12 and 13). 
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This is in agreement with earlier studies in the infinite temperature limit |31j . 

To summarize, xx and yy dynamic structure factors at low temperatures are not restricted to the two- 
fermion excitation continua and have (small) nonzero values outside these continua. They exhibit several 
washed-out excitations concentrated along the characteristic lines of the two-fermion excitation continua. 
The Dzyaloshinskii-Moriya interaction manifests itself in the constant frequency/ wave- vector scans influ- 
encing the detailed structure of such profiles. In the high-temperature limit xx and yy dynamic structure 
factors become K-independent. 

5 Conclusions 

In this paper, we have obtained the detailed dynamic structure factors S aa (K,u)), a = x, y,z of the spin-i 
anisotropic XY chain in a transverse field with the Dzyaloshinskii-Moriya interaction. The Dzyaloshinskii- 
Moriya interaction leads to nontrivial changes in the dynamic quantities. The two-fermion excitations which 
exclusively govern the zz dynamic structure factor, form three excitation continua and all of them manifest 
themselves even at zero temperature for sufficiently large strength of the Dzyaloshinskii-Moriya interaction. 
The two-fermion dynamic quantities have nonzero values in a restricted region of the k—u> plane; they may 
exhibit van Hove singularities (not only with exponent i but also with |); moreover, they may exhibit finite 
jumps at zero temperature. The xx and yy dynamic structure factors involve many-fermion excitations. 
However, the two-fermion excitations dominate their low-temperature behavior: at low temperatures these 
quantities show several washed-out excitation branches which correspond to specific lines of the two-fermion 
excitation continua. 

The Dzyaloshinskii-Moriya interaction clearly manifests itself in the frequency /wave- vector profiles which 
makes it possible to determine the magnitude of this interaction by measuring the dynamic structure factors. 
Dynamical structure factors can be measured by neutron scattering. Another experimental technique which 
yields dynamic quantities is electron spin resonance (ESR). If a static magnetic field along z axis and the 
electromagnetic wave polarized in a _L z direction are applied to the spin-i anisotropic XY chain with the 
Dzyaloshinskii-Moriya interaction, the experimentally measurable absorption intensity is given by 

/(wjocw S aa {0,LO). (5.1) 

Thus, the theoretical results for the dynamic structure factors presented in Section 4 are directly related to 
the ESR absorption intensity I(u>) (|5.1[) . Similar study for the spin-i isotropic XY chain has been recently 
reported in Ref. [J5] . A detailed analysis of the Dzyaloshinskii-Moriya interaction effect on ESR experiments 
on the materials, which can be modeled by the spin-^ anisotropic XY chain, seems to be an interesting issue. 
The work in this direction is in progress. 

Acknowledgments 

O. D. acknowledges the kind hospitality of the University of Magdeburg in the end of 2005 and in the 
beginning of 2006. T. V. acknowledges the kind hospitality of the University of Bayreuth in the spring and 
autumn of 2004. The paper was presented at the International Conference on Strongly Correlated Electron 



10 



Systems (Vienna, July 26th - 30th, 2005). O. D., T. V. and T. K. thank the Organizing Committee for the 
support for attending the conference. 

References 

[1] P. Lcmmcns, G. Giinthcrodt, and C. Gros, Phys. Rep. 375, 1 (2003). 

[2] M. Oshikawa and I. Affleck, Phys. Rev. Lett. 79, 2883 (1997); I. Affleck and M. Oshikawa, Phys. Rev. 
B 60, 1038 (1999). 

[3] J. Z. Zhao, X. Q. Wang, T. Xiang, Z. B. Su, and L. Yu, Phys. Rev. Lett. 90, 207204 (2003). 

[4] I. Tsukada, J. Takeya, T. Masuda, and K. Uchinokura, Phys. Rev. Lett. 87, 127203 (2001). 

[5] S. Bertaina, V. A. Pashchcnko, A. Stepanov, T. Masuda, and K. Uchinokura, Phys. Rev. Lett. 92, 
057203 (2004). 

[6] M. Kohgi, K. Iwasa, J. - M. Mignot, B. Fak, P. Gegenwart, M. Lang, A. Ochiai, H. Aoki, and T. Suzuki, 
Phys. Rev. Lett. 86, 2439 (2001). 

[7] P. Fulde, B. Schmidt, and P. Thalmeier, Europhys. Lett. 31, 323 (1995). 

[8] M. Oshikawa, K. Ueda, H. Aoki, A. Ochiai, and M. Kohgi, J. Phys. Soc. Jpn. 68, 3181 (1999). 
[9] H. Shiba, K. Ueda, and O. Sakai, J. Phys. Soc. Jpn. 69, 1493 (2000). 
[10] I. Yamada and H. Onda, Phys. Rev. B 49, 1048 (1994). 

[11] J. Pommer, V. Kataev, K. - Y. Choi, P. Lemmens, A. Ioncscu, Yu. Pashkcvich, A. Freimuth, and G. 
Giinthcrodt, Phys. Rev. B 67, 214410 (2003). 

[12] V. A. Ivanshin, V. Yushankhai, J. Sichelschmidt, D. V. Zakharov, E. E. Kaul, and C. Geibel, Phys. 
Rev. B 68, 064404 (2003). 

[13] H. - A. Krug von Nidda, L. E. Svistov, M. V. Eremin, R. M. Eremina, A. Loidl, V. Kataev, A. Validov, 
A. Prokofiev, and W. Afimus, Phys. Rev. B 65, 134445 (2002). 

[14] I. Yamada, M. Nishi, and J. Akimitsu, J. Phys.: Condens. Matter 8, 2625 (1996); V. N. Glazkov, A. I. 
Smirnov, O. A. Pctrcnko, D. M C K Paul, A. G. Vctkin, and R. M. Eremina, J. Phys.: Condens. Matter 
10, 7879 (1998); H. Nojiri, H. Ohta, S. Okubo, O. Fujita, J. Akimitsu, and M. Motokawa, J. Phys. Soc. 
Jpn. 68, 3417 (1999). 

[15] N. Shibata and K. Ueda, J. Phys. Soc. Jpn. 70, 3690 (2001). 

[16] T. Moriya, Phys. Rev. Lett. 4, 228 (1960); T. Moriya, Phys. Rev. 120, 91 (1960). 

[17] T. Antal, Z. Racz, and L. Sasvari, Phys. Rev. Lett. 78, 167 (1997). 

[18] T. Antal, Z. Racz, A. Rakos, and G. M. Schiitz, Phys. Rev. E 57, 5184 (1998). 

11 



[19] T. Antal, Z. Racz, A. Rakos, and G. M. Schutz, Phys. Rev. E 59, 4912 (1999). 

[20] M. D. Grynberg, T. J. Newman, and R. B. Stinchcombe, Phys. Rev. E 50, 957 (1994). 

[21] D. C. Mattis and M. L. Glasser, Rev. Mod. Phys. 70, 979 (1998). 

[22] D. N. Aristov and S. V. Maleyev, Phys. Rev. B 62, R751 (2000). 

[23] M. Bocquct, F. H. L. Esslcr, A. M. Tsvclik, and A. O. Gogolin, Phys. Rev. B 64, 094425 (2001). 

[24] M. Kenzelmann, R. Coldea, D. A. Tennant, D. Visscr, M. Hofmann, P. Smeibidl, and Z. Tylczynski, 
Phys. Rev. B 65, 144432 (2002); C. J. Mukherjee, R. Coldea, D. A. Tennant, M. Koza, M. Enderle, K. 
Habicht, P. Smeibidl, and Z. Tylczynski, J. Magn. Magn. Mater. 272-276, 920 (2004). 

[25] J. - S. Caux, F. H. L. Esslcr, and U. Low, Phys. Rev. B 68, 134431 (2003). 

[26] D. V. Dmitricv and V. Ya. Krivnov, Pis'ma v ZhETF 80, 349 (2004). 

[27] J. - S. Caux and J. M. Maillet, Phys. Rev. Lett. 95, 077201 (2005); J. - S. Caux, R. Hagemans, and J. 
M. Maillet, J. Stat. Mech.: Theor. Exp. P09003 (2005). 

[28] V. M. Kontorovich and V. M. Tsukcrnik, ZhETF 52, 1446 (1967) (in Russian). 

[29] Th. J. Siskcns, H. W. Capel, and K. J. F. Gaemers, Physica A 79, 259 (1975). 

[30] Th. J. Siskcns and H. W. Capel, Physica A 79, 296 (1975). 

[31] J. H. H. Perk and H. W. Capel, Physica A 92, 163 (1978). 

[32] O. Derzhko and A. Moina, Ferroelectrics 153, 49 (1994); O. V. Derzhko and A. Ph. Moina, Condensed 
Matter Physics (L'viv) 3, 3 (1994). 

[33] A. S. T. Pires, J. Magn. Magn. Mater. 223, 304 (2001). 

[34] O. Derzhko and T. Verkholyak, Proceedings of Institute of Mathematics of NAS of Ukraine 50, 692 
(2004); O. Derzhko and T. Verkholyak, Czech. J. Phys. 54, D531 (2004); O. Derzhko and T. Verkholyak, 
Physica B 359-361, 1403 (2005). 

[35] The difference between the one-dimensional isotropic and anisotropic XY models was also underlined 
in: J. H. H. Perk and H. W. Capel, Phys. Lett. A 58, 115 (1976). 

[36] O. Derzhko and T. Krokhmalskii, Phys. Rev. B 56, 11659 (1997); O. Derzhko and T. Krokhmalskii, 
physica status solidi (b) 208, 221 (1998); O. Derzhko, T. Krokhmalskii, and J. Stolzc, J. Phys. A 33, 
3063 (2000); O. Derzhko, T. Krokhmalskii, and J. Stolzc, J. Phys. A 35, 3573 (2002). 

[37] D. Karevski, Relaxation in Quantum Spin Chains: Free Fcrmionic Models, in Order, Disorder and 
Criticality. Advanced Problems of Phase Transition Theory, edited by Yu. Holovatch (World Scientific, 
Singapore, 2004), pp.67-107. 

[38] J. H. Taylor and G. Miiller, Physica A 130, 1 (1985). 



12 



[39] H. Suzuura, H. Yasuhara, A. Furusaki, N. Nagaosa, and Y. Tokura, Phys. Rev. Lett. 76, 2579 (1996). 
[40] Yongmin Yu, G. Miiller, and V. S. Viswanath, Phys. Rev. B 54, 9242 (1996). 

[41] O. Dcrzhko, T. Krokhmalskii, J. Stolze, and G. Miiller, Phys. Rev. B 71, 104432 (2005); O. Derzhko, 
T. Krokhmalskii, J. Stolze, and G. Miiller, Physica B 378-380C, 445 (2006). 

[42] O. Derzhko, T. Verkholyak, T. Krokhmalskii, and H. Biittner, Physica B 378-380C, 443 (2006). 

[43] G. Miiller, H. Thomas, H. Beck, and J. C. Bonner, Phys. Rev. B 24, 1429 (1981). 

[44] I. Affleck, Field theory methods and quantum critical phenomena, in Fields, Strings and Critical Phe- 
nomena, edited by E. Brezin and J. Zinn-Justin (Elsevier, Amsterdam, 1989), pp. 563-640; J. von Delft 
and H. Schoeller, Annalen der Physik 7, 225 (1998). 

[45] Y. Maeda and M. Oshikawa, Phys. Rev. B 67, 224424 (2003). 



FIGURE CAPTIONS 

FIG. 1. Gray-scale plot for the zz dynamic structure factor S zz (k,ui) (|2.3| for the spin chain (jl.ip with 
J = 1, 7 = 0.5 at low temperature, (3 = 50, for different strengths of the Dzyaloshinskii-Moriya interaction, 
D = (right column), D = 0.5 (middle column), D = 1 (left column), and different fields, O = 0, 0.5, 1, 
(from bottom to top). 

FIG. 2. Frequency profiles of S zz (k,uj) at different wave- vectors k = (left panel), k — -| (middle panel) 
and k = 7r (right panel) for the spin chain (|1.1[) with J = 1, 7 = 0.5, D — (dotted lines), D = 0.5 (dashed 
lines), D = 1 (solid lines), SI = 0.5 at low temperature (3 = 50. 

FIG. 3. The zz dynamic structure factor (|2.3[) for the spin chain (|1.1[) with J = 1, 7 = 0.5, D = 1, 
= 0.5 at different temperatures (3 — 10 (a), /3 = 1 (b) and (3 = 0.1 (c). 

FIG. 4. The gray-scale plots of S zz \hl,uj) (a), S z 2 z (n, to) (b) and S z Z z \k,lo) (c) for the set of parameters 
J = 1, 7 = 0.5, D =1,Q. = 0.5, (3 = 50. 

FIG. 5. For comparison s'i) (k,lo) (|3.ip with B^\ki,k) = 1 for the same set of parameters as in Fig. 4. 
a: j = 1, b: j = 2, c: j = 3. 

FIG. 6. Towards the properties of the two-fermion excitation continua; J = 1, 7 = 0.5, D = 1, Q, = 0.5; 
left panels - infinite temperature limit T — > 00 ((3 = 0), right panels - zero temperature limit T = ((3 — > 00); 
j = 1 (panels a and d), j = 2 (panels b and e), j = 3 (panels c and f). 

FIG. 7. Frequency profiles (with 5^ («!,«) = 1) for the chain COD with J = 1, 7 = 0.5, D = 1, 

il = 0.5, (3 = which illustrate van Hove singularities with different exponents | and |. Panel a: S^(k,uj) 
(13. 8|) vs tj at k = 0.5 (solid line), k = 2 (dashed line), n = 3 (dotted line); panel b: S'^^w) (|3.8|) vs a; at 
k = 0.9 (solid line), k 1.07844531 (dashed line; the low-frequency singularity has power-law exponent |), 
k = 1.2 (dotted line). 
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FIG. 8. The xx dynamic structure factor S xx (k, lu) for the spin chain with J = 1,7 = 0.5 at low tem- 
perature /3 = 50. D = (right column), D — 0.5 (middle column), D = 1 (left column), f2 = 0, 0.25, 0.5, 1 
(from bottom to top). 

FIG. 9. The yy dynamic structure factor S yy {n, co) for the spin chain (jl.ip with J = 1,7 = 0.5 at low tem- 
perature (3 = 50. D — (right column), D = 0.5 (middle column), D — 1 (left column), f2 = 0, 0.25, 0.5, 1 
(from bottom to top). 

FIG. 10. Frequency profiles of S xx (k,uj) at k — (left panel), k = § (middle panel) and k — n (right 
panel) for the spin chain (jl.ip with J = 1, 7 = 0.5, Z? = (dotted lines), D = 0.5 (dashed lines), D = 1 
(solid lines), f2 = 0.5 at low temperature (3 = 50. 

FIG. 11. Frequency profiles of Syy(K,uj) at k — (left panel), k = ^ (middle panel) and k — 7r (right 
panel) for the spin chain (jl.ip with J = 1, 7 = 0.5, D = (dotted lines), I? = 0.5 (dashed lines), D = 1 
(solid lines), = 0.5 at low temperature f3 = 50. 

FIG. 12. The xx dynamic structure factor S xx (n, to) for the spin chain (jl.ip with J = 1, 7 = 0.5, D = 1, 
fl = 0.5 at different temperatures: (3 — 10 (left panel), /3 = 1 (middle panel) and (3 — 0.1 (right panel). 

FIG. 13. The yy dynamic structure factor S yy (K,io) for the spin chain (|l.ip with J = 1, 7 = 0.5, D = 1, 
il = 0.5 at different temperatures: (3 = 10 (left panel), /? = 1 (middle panel) and (3 = 0.1 (right panel). 

FIG. 14. S' 22 (K,ti;), S xx (k,u>) and S yy (K^uj) (from top to bottom) for the spin chain (jl.ip with J = 1, 
7 = 0.5, Z) = 1. O = 0.5 at low temperature (3 = 50. 
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